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. abstract 

a.: 

, In an accompanying paper, we have formulated two types of regularization methods to calculate the scalar self-force on a 

particle of charge q moving around a black hole of mass M jl), one of which is called the "power expansion regularization". 
In this paper, we analytically evaluate the self-force (which we also call the reaction force) to the third post-Newtonian (3PN) 
order on the scalar particle in circular orbit around a Schwarzschild black hole by using the power expansion regularization. It 
is found that the r-component of the self-force arises at the 3PN order, whereas the t- and (^-components, which are due to the 
radiation reaction, appear at the 2PN and 1.5PN orders, respectively. 
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Thanks to the recent technological advance, we have almost come to the stage that gravitational waves are de- 
tectable. The observation of gravitational waves is absolutely a new window to observe our universe. We also expect 
that the observation of gravitational waves provides a direct experimental test of general relativity. 

There are several on-going projects of gravitational wave detection in the world, such as LIGO [Q], VIRGO [[|, 
GEO Q and TAMA H. In addition, LISA M is proposed a future project. LISA is unique in that it is a space-based 
^Jy observation. Since the detector will be launched in space, it will be free from the seismic noise, and it will be remarkably 
sensitive to the low frequency gravitational waves below 1 Hz. In contrast to the ground-based gravitational wave 
detectors whose main targets are inspiralling binaries of solar-mass compact objects (either black holes or neutron 
stars), LISA will detect the gravitational waves from a solar- mass compact object orbiting a supermassive black 
hole at galactic centers. In the latter situation, the radiation reaction effect over a few orbital periods will be much 
smaller than the case of compact binaries because of the large mass ratio between the black hole and the compact 
object. However, the accumulated "secular" effect of radiation reaction will be as important as in the case of compact 
binaries. Furthermore, the spacetime around a supermassive black hole is expected to be described by the Kerr 
geometry with non-negligible spin. Hence the relativistic effect, especially the effect of the non-spherical geometry, 
will significantly affect the orbital evolution. In this situation, the calculation of the orbital evolution by means of 
the standard post-Newtonian expansion method seems formidable to perform. Instead, the black hole perturbation 
approach seems much more relevant, where the compact object is approximated by a particle orbiting a Kerr black 
hole (and perturbing the geometry). 

For a binary system of almost equal mass objects, the orbit is expected to become sufficiently circular at the final 
inspiral stage (~ 10 5 cycles before coalescence) by the radiation reaction. However, this will not be the case for a 
binary systems of extreme mass ratio, and the orbit may not become circular until the very last moment. Thus we 
are required to consider general orbits. 

Previous work on radiation reaction in the black hole perturbation approach has been based on conservation of 
the total energy and angular momentum of the system (see Ref. Q| for a review). However, this approach cannot 
work for general orbits around a Kerr black hole, since the orbit is characterized by the three parameters; the energy, 
angular momentum (with respect to the symmetry axis), and the Carter constant Since the Carter constant is not 
associated with the Killing vector field of the Kerr geometry, its evolution cannot be determined by the conservation 
laws. Instead, it is necessary to derive the reaction force acting locally on the orbit. 



Recently, the equation of motion of a point particle on a curved background with the effect of gravitational radiation 
reaction was derived by Mino, Sasaki and Tanaka by solving the Einstein equations by means of matched asymptotic 
expansion ||, and by Quinn and Wald by an axiomatic approach p0| . 

The essential problem of the equation of motion of a point particle was that the metric perturbation induced by the 
particle diverges at the location of the particle. Apparently, divergence of the metric perturbation means the break 
down of the linear approximation. Nevertheless, we may expect that the strong equivalence principle still holds to 
an extent, and that the notion of the equation of motion of a point particle is still useful. The work of Ref. Q and 
Ref. jlO) has shown this is indeed the case, and that the extra force by particle's self-gravity can be derived from the 
linear metric perturbation by subtracting the singular part appropriately. 

The singular part irrelevant to the motion can be exactly identified by the covariant Hadamard method Jll]] as the 
coincidence limit of the part that directly propagates to field points along the background null cone of the source. 
The rest of the metric perturbation is called the tail part, which propagates within the null cone due to the curvature 
scattering. In practice, the direct part can be calculated by using the method of DeWitt and Brehme On the 
other hand, the calculation of the tail part requires full knowledge of the metric perturbation, and it is impossible to 
calculate it for a general spacetime. Fortunately, in the case of the black hole perturbation, it is in principle possible 
to calculate the full metric perturbation |l3) by using the Regge- Wheeler- Teukolsky formalism [Q. Then we may 
extract the tail part by subtracting the direct part from the full metric perturbation. 

Here, however, we must mention a couple of non-trivial obstacles when we attempt to perform the extraction 
procedure: 

1) Subtraction problem 

In the Regge- Wheeler- Teukolsky formalism, the full metric perturbation is calculated in the form of Fourier- 
harmonic series.^] This fact totally obscures the local spacetime behavior of the metric perturbation around the 
particle. On the other hand, the direct part is given only locally around the location of the particle. This nature 
of locality makes it very difficult for us to transform the direct part to the Fourier-harmonic series form. Thus 
the subtraction of the direct part from the full metric perturbation is a highly non-trivial task. 

2) Gauge problem 

The Hadamard prescription which identifies the direct part is, by definition, applicable only to hyperbolic 
differential equations, and specifically to the perturbation equations in the Lorentz (harmonic) gauge in the case 
of our interest. On the other hand, the metric perturbation directly obtainable in the Regge- Wheeler- Teukolsky 
formalism is not the one in the Lorentz gauge. Hence an appropriate gauge transformation should be worked 
out. 

In addition to the above, in the case of a Kerr background, which is of our ultimate interest, there is a much more 
intricate technical issue of how to treat the spheroidal harmonics: Neither the eigenfunction nor the eigenvalue has 
simple analytical expressions, and they are entangled with the frequency eigenvalues of Fourier modes. Although our 
ultimate goal is to overcome these difficulties altogether, since each one is sufficiently involved, we choose to proceed 
step by step, and tackle the subtraction problem first. 

As a solution to the subtraction problem, we have recently formulated two types of regularization methods [0, 



extending the earlier work [[15 16 1. We termed the two methods the 'power expansion regularization' and the 'mode- 
by-mode regularization'. Although the methods we have developed should be applicable to the gravitational case in 
principle, we have focused on the scalar case since it is free from the gauge problem: The Regge- Wheeler or Teukolsky 
equation reduces to the simple scalar d'Alembertian equation. As for the mode-by-mode regularization, a variant of it 
adapted to the gravitational case has been applied to a particle in circular orbit around a Schwarzschild black hole in 
Ref. JlTj] successfully but only to the 1PN order. The purpose of the present paper is to demonstrate the effectiveness 
of the power expansion regularization explicitly and analytically by applying it to the simplest case of a scalar particle 
in circular orbit around a Schwarzschild black hole to the 3PN order. 



Analytic solutions of the Teukolsky equation were derived by Mano, Suzuki and Takasugi |L8 19 in the series of 
hypergeometric functions. We use their result to construct the Green function to the 3PN order where the technical 
advantage of the power expansion regularization becomes clear. The direct part has been obtained under the local 
expansion for general orbits in Ref. Ill . 

There have appeared some papers discussing the extraction of the scalar self- force pp|-^5| . Recently, using the mode- 
sum regularization prescription proposed in Ref. p3], the self- force acting on an electric or scalar charge has been 



"Throughout the paper, the Fourier modes refer to e and the harmonic modes to Ye m (8, <f>). 
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discussed in the spacetime of spherical shells. Recently the nature of the gravitational self-force has been discussed 
by Detweiler |p6| . 

This paper is organized as follows. In Sec. 0, we briefly review the power expansion regularization developed in 
Rcf. Q. In Sec. Ill, we consider the full Green function to the 3PN order in the Fourier-harmonic expanded form, and 



perform the power expansion. In Sec. [Tv| , we perform the power expansion of the direct part evaluated by the local 
expansion. Then in Sec. 0, we derive the regularized scalar self-force by applying the power expansion regularization 
to the results obtained in the previous two sections. Finally, we summarize our result and discuss possible extensions 



of our work in Sec. Nil 



II. METHOD OF REGULARIZATION 

The reaction force on a point scalar charge is given by 

F°(j)= lim (F a [<f^ u }(x) - F a [<f> dil ](x)) = lim F"[0 tail ](x) , (2.1) 

x — >z{t) x — >z(t) 

where (j) lul1 is the scalar field induced by the point source, dlr is the direct part of the field defined in Ref. ]2^ , p^ , [l0| ,p| , 
^taii _ ^fuii _ 0dn^ an( j jr"[___] i s a linear tensor derivative operator to derive the force from the scalar field. To obtain 
</> fu11 , we construct the full Green function in the form of Fourier-harmonic series. Since both the full scalar field 
and the direct part diverge at the location of the particle, we calculate them at a field point off the particle where 
the subtraction of the divergent part is well-defined. We call this field point xf eg the regularization point. There is 
freedom in the choice of the regularization point. Let the location of the particle at its proper time r = tq be 

{z Q } = {z a (r )} = {t , r , 6 , <M • (2.2) 

For the power expansion regularization, which will be briefly explained later in this section, we take {x" og } = 
{t,r, 6> ,<M- 

We consider a point scalar charge q moving in the Schwarzschild background, 

ds 2 = - (l - ™\ dt 2 + (l - 1 dr 2 + r 2 (d9 2 + sin 2 9d<p 2 ) , (2.3) 

where {x a } = {t, r, 9, 4>} are the Schwarzschild coordinates, and M is the black hole mass. The charge density of the 
scalar particle in circular orbit is given by 

^..L jV^W) , (2 .4) 



where z a is the trajectory of the particle, 



z -(r) = («V, r , -, v*t\ ; v* = J r ° ", W * = - A /^^^, (2.5) 

and v a is the four velocity of the particle. We have assumed the orbit is on the equatorial plane without loss of 
generality. 

The full scalar field is calculated by using the Green function method as 

full (x) = qj drG full (x, z{t)) , (2.6) 
where the Green function satisfies the wave equation 

V<*V a G™(x,x>)= S{4){ ^Z_f\ (2.7) 

and we impose the retarded boundary condition on G full (x,x'). The full Green function is decomposed into the 
Fourier-harmonic modes as 
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G™(x,x') = f —e^-^^MY^m^^) 



(2.8) 



Then the wave equation (2.7) reduces to an ordinary differential equation for the radial Green function as 



2M\ d 2 2(r - M) d 
r I dr 2 r 2 



dr 



£(£+l) 



V 



2M 



(2.9) 



We construct the radial function of the full Green function by using homogeneous solutions of ( |2.9[ ) which can be 
obtained by a systematic analytic method developed in Ref . |Q . 

On the other hand, the local expansion of the full Green function gives the direct part of the scalar field, 



</> dir (z) =Q / dTG d "(x lZ { T )) 



where G is the direct part of the Green function. It is given in a covariant manner as 

1 



G dir (x,x') = -— 6>[£(x), x'} y/A(x, x')<%(x, x 1 )) 



(2.10) 



(2.11) 



where <t(x,x') is the bi-scalar of half the squared geodesic distance, A(x,x') is the generalized van Vleck-Morette 
determinant bi-scalar, X(x) is an arbitrary spacelike hypersurface containing x, and 0[£(x),x'] — 1 — (9[x',£(x)] is 
equal to 1 when x' lies in the past of £(x) and vanishes when x' lie in the future. 

We now briefly descr ibe the power expansion regularization procedure. As mentioned earlier, for a given point on 
the particle trajectory (2.2), we choose the regularization point as 



{< e J = {t,r,0 o ,<fa} ■ 



(2.12) 



which is assumed to be on a suitably chosen spacelike hypersurface: See Eq. (3.8) in Ref. Q. It should be noted that 
t is a function of r. In the following, we take tq = for simplicity. 

Assuming (r — ro)/ro <C 1, we expand the full scalar field and its direct part in powers of the radius r as 



jfull 



(x rcg ) = rV ull(n) , dir (x rcg ) = ]T rV 



dir(n) 



(2.13) 



Similarly, for the calculation of the self-force, we expand V p 



/.full 



and V, 



as 



V^ full (x reg ) = r n ^ U{n) , V^ dir (x rcg ) = ^"< ir 



(n) 



(2.14) 



For 



/.full 



and V^'" 11 , this expansion is done first for each spherical harmonic mode, followed by summation of the 
Fourier- harmonic components to obtain full (™) and cj/^ 1 ^ n The heart of the power expansion regularization resides 



in the fact that the coefficients </> full (™) and 



,full(n) 



as well as 6 A " {n) and 



i dir(n) 



of the power series are finite, though 



these series as a whole diverge in the limit r — > tq. Hence we may extract the tail part as 

tail (x rcg ) = r"</) tail(n) ; tail(Tl) = full(Tl) - (t> diT{n) . 



V^ tail (x rcg ) = r n ^W ; < ail(ri) = < U(n) - C (n) ■ 



(2.15) 



The tail part of the scalar field is obtained by summing over n and taking the limit r — > ro . This extraction procedure 
is called the power expansion regularization. 



III. POWER EXPANSION OF THE FULL SCALAR FIELD 



We first consider the full Green function decomposed in the Fourier-harmonic modes (2.8), The retarded radial 
Green function becomes 



4 



^ (r ' r,) = liU^) (^( r )#™(*W - r) + WCr - r')) ; 



where 0"^^ is a homogeneous solution which vanishes on the past horizon (when multiplied by e~ lu)t ), and 4>^ luJ 
is a homogeneous solution which vanishes at the past null infinity |2§f| . They are called the in-going and up-going 
solutions, respectively. 



Homogeneous solutions of (2.9) were fully investigated in Ref. |lq,tt9[ and we have various analytic expressions for 



the homogeneous functions. In this paper, we consider the slow-motion expansion, in which we assume 

wr k d, ujM k v 3 ; v < 1. (3.2) 

Detailed properties of the radial Green function for the calculation up to the 3PN order are described in Appendix 
[A|. For technical reasons, we classify the harmonic components of the Green function into the three parts, 

G$g>(x,x')= f gijU-M*^ J2 9t™\ry)Y im (e,m m (0', </>'), (3-3) 

G^ ) (x,x')= f ^iV-C*-0 dti\r,r>)Y em (6^)Y; m (e>^>), (3.4) 

J n £>2,rn 

G^ Q<1) {x,x') = /^O e " Ml "° E &(rS)Y tm (0,<l>)Y£ n (0',<t>'). (3.5) 

J ^ 2=0,1, rn 

The suffices (sym) and (rad) represent that they are the time-symmetric and radiative parts, respectively, of the Green 
function. That is, we write the (retarded) Green function as 

^»full ^f(sym) _j_ G (rad) . 

G (sym) = _( G (ret) + G (adv)) > G (rad) _ ^(ret) _ G (adv)j ; (3 g) 

where G^ adv ^ stands for the advanced Green function. Correspondingly, the full scalar field is divided into the three 
parts as 

<fe } 0=) = / d±x'^\W)Gf^Ax,x')p{x') , (3.7) 



«>2) W - I «* * V lyiv^ /^(£>2) 



^> 2 \W = / dV^^G^Or, £>(*'), (3.8) 



<$=o,i)(*) = rfVV^G^^^')^')- (3.9) 

The division of the Green function to the harmonic modes of £ > 2 and £ = 0, 1 is purely technical because of some 
non-systematic analytical behavior of the modes £ — 0, 1 as discussed in Appendix |X| On the other hand, there is a 
physical reason for dividing the Green function to the time-symmetric part and the radiative part. It is known that 
the radiation reaction to the energy and angular momentum of the particle (i.e., any conserved quantity associated 
with a background Killing vector field) can be calculated by using the radiative Green function G^ rad ^ [|{| . In our case 
of circular orbit around a Schwarzschild black hole, we therefore expect no contribution from the time-symmetric part 
to the t- and (^-components of the reaction force. Nevertheless, one cannot simply discard the time-symmetric part 
since the tail part is not identical to the radiative part. In particular, the r-component of the reaction force, which 
does not contribute to the energy and angular momentum loss in our case, cannot be obtained from the radiative part. 
Since our purpose here is to demonstrate the validity of the power expansion regularization, we regularize G^>™j to 
show the vanishing of the t- and (^-components of the reaction force due to the time-symmetric part explicitly, as well 
as to derive the non- vanishing r-component of the reaction force. As for G^^ and G^ ^, they are always finite. 
Hence the power expansion is unnecessary for these parts. 

The time-symmetric part of the radial Green function is expanded in powers of r as 

&\r,r>) = £ (6(r> r)^^(r') + 9(r - rO^C^V)) . (3.10) 
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We note that because of the singularity at r = ro, we have two different forms of the power expansion at < r < tq 
and ro < r. The final regularized force, however, is the same for either choice of the regularization point. 

For the calculation to the 3PN order, it is found that we only need g^J for — 3 < n < 6 and for — 6 < n < 3. 

The coefficients g™^ and g^^ are summarized in Appendix [b|. After summing over u>, £ and m, we obtain the 
power expansion of the time-symmetric Green function as 

r' n i \ \ 

+e(r-r')— I G^(0^y,6 , ,<l> , ;u 1 )), (3.11) 



where 



G^ n) (9A;r',6',<f>';u,)= £ ^^"V^M)^, 0') . (3.12) 

«>2 



One sees that only a finite number of harmonic modes contribute to each expansion coefficient Q ln / out ( n ) because of 
the truncation of the terms at large |n — £\, which grea tly eases the calculation. 
Now the time-symmetric part of the scalar field (|3.7D becomes 



= f E ( - r )^ Mn) fa) + °(r ^^^(vo)) , (3.13) 

where 

0i n/out(„) =^g^ n - e \ ro )Y lrn (6,ct>)Y; m (n/2,0)e- imm , (3.14) 

Cm 

and we have introduced the angular frequency fl — v^/v*. Taking the derivative of the scalar field, we have the power 
expanded scalar field at the regularization point as 

Vt^j^eg) =0, (3.15) 

Vr4$$M = 4E (V° ~ r)n^0fW(n,) - 9(r - r )(n + , (3.16) 

V^g(z reg )=0, (3.17) 



V^^> 2 ) (*»«) = -^M^j (*«*) = , (3.18) 



where 



c( out(n) = E^ / :: t ^ ) (^)ir, m (-/2,o)i 2 . (3.i9) 

tm 

The fact that the ^-component vanishes comes from our assumption of circular orbit, while the reason for the vanishing 

t- and (^-components is a bit subtle. It comes both from t he a ssumption of circular orbit and the symmetry of 
under u <-» — lu. The formula used for the m-sum in Eq. (3.19) is presented in Appendix ^|. 



IV. POWER EXPANSION OF THE DIRECT PART 



Here we evaluate the derivative of the direct part of the scalar field under the local expansion. The details of this 
evaluation procedure are discussed in Ref. Q. By calculating Eqs. (3.15-7) in Ref. Q, we obtain the regularized direct 
part of the scalar reaction force as 



G 



V t dir (x rcg ) = O(r-r o ), (4.1) 
di q f 2M sgn(r-ro) / 1 M 2 (r -M) , 

VA M ^^y 1 - — ~v^w V ~ 4 Wo - mwo - m {r - ro) 

+0(r - r ) , (4.2) 

V e dir (x rcg ) = O(r-r o ), (4.3) 

(x rcg ) = 0(r - r ) , (4.4) 



dir 



where sgn(r — ro) = ±1 for r — ro ^ 0, and we have ignored the terms of order r — rg as they do not contribute 
to the force in the coincidence limit r — > ro. It should be mentioned that the above are exact in the sense that no 
slow-motion approximation is employed. Note that the second term proportional to M /Vq in the parentheses on the 
right-hand-side of the radial component gives a manifestly finite contribution in the limit r — > ro. Of course, these 
simple forms of the direct force components owe to the assumption of circular orbit: They will be more complicated 
if the radial velocity is non-zero. 
Using the formula, 

^ = £ r„)(n + 1)-^ + 9(r - r)(n + 1) — , (4-5) 

(r-r ) 2 f-^ Q \ r»+ 2 r J + V 

we finally obtain the power expansion of the direct part as 

V^ dir (z reg ) = 0, (4.6) 



V^ dlr (a; reg ) 




M 2 (r -M) 



4r 2 (r -2M) 2 (r - 3M) 
1 M 2 (r -M) 



(4.7) 



'rj+ 2 4r 2 (r - 2M) 2 (r - 3M) 

V e dir (^ O g) -0, (4.8) 
V^ dir (z reg ) = . (4.9) 

We find that only the r-component of the force diverges in the limit r — > ro. 

V. EXTRACTION OF THE SCALAR REACTION FORCE 

Now we are ready to discuss the extraction of the tail part under the power expansion regularization. Since the 
r-component of the full (bare) force is the only one that diverges in the coincidence limit, we first consider the simpler 
cases of the t-, 9- and (^-components, and separately treat the r-component afterwards. 

For the ^-component, all the parts of the force vanish identically because of the assumption of circular orbit, 

V^ tail Oz) = 0. (5.1) 
For the t- and (^-components, which are equal to each other mod ulo th e fact or — 1/17, both the time-symmetric part and 



the direct part of the force vanish identically, as given in Eqs. ( 3.15 ) and ([h9j). Hence the only possible contributions 
come from Vt^fe^ and Vt0(^ = o,i), in accordance with our expectation. We find the former contribution arises at 
the 3PN order and the latter at the 2PN order, which are added together to give 

v t tail (z) = Vt^gj(«) + Vt^,^*) 

q (\ 4 46073 fi _ 7 A . „, 

' v 4 + — — v 6 + 0{v 7 ) , (5.2) 



47rr 2 \3~ ' 30 
V^ tail (z) = -I Vt tail (z) 

q fl , 46073 



-v 



47rro V3 30 



v 5 + 0{v 7 )) . (5.3) 
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It is easy to check that the leading terms of the above are consistent with that of the energy loss rate by scalar 
radiation evaluated from the energy-momentum tensor of the scalar field at future null-infinity. The leading terms 
are also in agreement with Gal'tsov p9| . 

Now we turn to the r-component of the reaction force. Sinc e on ly the subtraction of the divergent terms is non- 
trivial, taking account of the form of V r <f> dir (x ieg ) given by Eq. (4.2), we introduce the auxiliary coefficients f 1D / ont ( n ) 

in the extraction procedure of the tail part from V r 0^>™j (^rcg) as 

2M sgn(r — 7*0) 



\t>2) 



( X rcg) 



1 



To 



(r - r ) 2 
r.n-1 
7+1 



Then we obtain 



f in(l 
/ in(2 

f n(3 
yin(4 

yin(5 

yin(6 

f in(7 

yin(n 

y*out(n 
y*out(— 4 

yout( — 3 

yout(-2 
yout( — 1 
yout(0 

y*out(l 

yout(2 

yout(3 

yout(4 

y*out(n 



0(« 7 ) 
1 + « 2 

33 4 

w 

20 



(n < 1) , 

ii v 4 + 331^6 + Q ,7y 

10 100 v ; 

173713 



14700 



v 6 + 0(v 7 ). 



3 2 ^ 4 
-1? 1> 



10" 10 
30154441 



215419 
22050 



v 6 + 0{v 7 ) , 



9604980 



v 6 + 0(v 7 ) , 



1 4 237035 

56 V + 30918888 

44372 6 _ 7 . 

V + 0(v 7 ) , 



w 6 + 0(v 7 ) , 



10735725 
1288333 



0(v 7 



527482800 

-6n(939 - 912n - 224n 3 + 224n 4 - 1972n 2 ) 
(2n - 7)(2n - 5)(2n - 3)(2ra - l) 2 (2n + 1) 2 (2 + 3) 2 '' 
0(v 7 ) (n < -4) , 
v 6 + 0(t; 7 ) , 

~v 6 + 0(v 7 ), 
15 

-^ 4 + ^ 6 + <V), 
+<V), 

3 



+ 0(v 7 ) (n>7) 



7 



7079 



1+ 2^-8^- 5040 u6 + ^' 
5 35 69300 v ; 



9 58 .. 
-6?; 2 + — w 4 + 



3686944 



525525 



+ 0(w 7 ) , 



12656296 



v 6 + 0{v 7 ) 



where v = flro (and flM = r 



72 4 ^ 
U + 

5 429975 
3396569326 fi ^, 7 . 

106135029 v ; 

-6(n + 1)(327 + 1120n 3 + 224n 4 - 1464n + 44n 2 ) 
(2n - l) 2 (2n + l) 2 (2n + 3) 2 (2n + 5)(2n + 7){2n + 9) 



+ 0(v 7 ) (n>4) 



(5.4) 

(5.5 
(5.6 

(5.7 

(5.8 

(5.9 

5.10 

5.11 

5.12 

5.13 

5.14 
5.15 

5.16 

5.17 
5.18 
5.19 

5.20 

5.21 

5.22 

5.23 

5.24 



As for the parts V r < 



,(rad) 
\l>2) 



and V r <; 



(^=1,0) ! 



we find the former is higher than the 3PN order and only the latter 
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contributes at the 3PN order. In fact, the radiative Green function does not contribute to the radial force for a 
circular orbit to full ord er (see Appendix |a|) . 

By summing up ( |5.4| ) and adding the finite terms in the direct part together with V r 0(£ = i j o), we finally find 

9 ^ u 6_^ 7r V + -u 6 ln(u)+0( W 7 )V (5.25) 



Anvrl V 12 64 3 

As expected, both choices of the regularization point, either at r < ro or r > ro, give the same reaction force in the 
limit r — » r . It is noted that the above result is consistent with that of the numerical calculation in Rcf. pl|. 



VI. CONCLUSION 

In a separate paper |l| , we have proposed two technical methods to derive the scalar reaction force on a particle in 
the Schwarzschild background, the power expansion regularization and the mode-by-mode regularization. We expect 
these methods will be useful also for the case of gravitational reaction force. In fact, a variant of the mode-by-mode 
regularization method has been successfully applied to the gravitational case, although only the 1PN calculation for 
circular orbit has been done (lTj . 

In this paper, we have demonstrated the usefulness of the other regularization method, the power expansion 
regularization method, by applying it to a scalar particle in circular orbit around a Schwarzschild black hole and 
calculating the scalar reaction force to the 3PN order. In the power expansion regularization method, the tail part 
of the scalar field, which is responsible for the reaction force, is obtained by subtracting the direct part from the full 
scalar field at a regularization point specified by the radial coordinate r, followed by taking the limit r — > ro where ro 
is the radial coordinate of the particle. We have found that the r-component of the scalar reaction force arises at the 
3PN order, and the t- and ^-components at the 2PN and 1.5PN orders, respectively. 

Since the derivation of reaction force for non-circular orbits is very important for binary systems of extreme mass 
ratio, and the gravitational waves from such binaries are expected to be detected by future space-based gravitational 
wave detectors such as LISA, it is indispensable for us to proceed to the case of general orbit. In this respect, although 
the calculation is bound to become more complicated when we consider non-circular orbits where the radial motion 
is present, the power expansion regularization method seems to be still applicable. A general framework for general 
orbit in terms of the power expansion regularization scheme has been developed in Ref . jlj . As a next step, we plan 
to consider the eccentric or the radial orbit with these regularization methods. 
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APPENDIX A: FULL GREEN FUNCTION 



In this appendix, we summarize the derivation of the full radial Green function. Here the in-going and up-going 
homogeneous solutions are denoted by 0f n and <f)^ p , respectively, where v is called the renormalized angular momentum 
Q,^9[, and is equal to t in the limit wM — -> 0. 

The radial Green function is expressed as 

r ') = w L ^ \ (CMfllpCfW - r) + « p (r)CMJ(r - /)) ; 

WWCCp) = (l - ^) (^cm)^pW - (J:#p(r))c(r)) • (Ai) 







For technical convenience, we express the homogeneous solutions 4>^ n and c/£ p in terms of the solutions <f>" and 4> c v 1 , 
which are expressed in terms of a series of the Coulomb wave functions |18| , as 

0Kp=7^c+Mc"~ 1 . (A2) 
where the properties and the relations of the coefficients {a„, V , j v , 8 U } have been discussed extensively in 
Ref. Hill. The function <p% is denoted by R v c in Ref. |||. It may be noted that 4>%(r) oc j t (wr) in the limit 
e — » 0, where ji is the spherical Bessel function, and we may choose its phase so that <f>" is real. 
We then divide the Green function to the time-symmetric and radiative parts, 

o fu11 (r r'\ — « (sym) lV r'\ A- J 13 ^ (r r'\ 
9lmui l r : r ) — 9hnu) V ' r ) + 9tmw V ' T I ' 

where 



(A3) 



i 



-5R 



1 - A/> 

^(r)^(r') + 5R 



. 1 - Pv% 

+9(r — r')\u «-> — v — 1 , /3„ <-> 7^ 



1 - 



(A4) 



1 - $y% 



1 - /9i/7i> 



C(r)C(r')+9 



1 - PvJv 



(A5) 



In the above, we have assumed av ^ and 8 V ^ 0, and introduced the coefficients {/3 V , j u } := {P v /a y , -fu/Su}- Using 
the result obtained in Ref. Jig ], we find the coefficients {/3„, 7^,} behave under post- Newtonian expansion as 

R|)9„] = 0(« 6£ + 3 ) , = 0(« M ) , »[>] = 0(« 6 ) , = (-l) f+1 . (A6) 

Although the analysis in Ref. [ p~8|Jl9| is very useful, it turns out that the post Newtonian expansion of <^ v ~ x breaks 
down at the (£ + 2)th PN order. Thus for our analysis to the 3PN order, the modes I = 0, 1 need a separate treatment. 

So, we consider the three cases, gf^^ and g'f^j for £ > 2, and gi muJ for I = 0, 1, separately. 



1. Time-symmetric part for £ > 2 

To the 3PN order, inspection of the post-Newtonian orders of {/3 V , y„} shows that the time-symmetric Green 
function for i > 2 can be approximated as 

&\rS) = w ( l ,-,-1, W rW c {r)r c v ^{r') + 8(r r')^ (rW c (r')} + ■ • • ■ (A7) 
Wtmw {((>%, <Pc ) 

The homogeneous solution is expanded to 0(u 6 ) as 

^ z 2 It z 4 {£ 2 -M-\0)ez £{£-l) 2 e 2 



<K{z) = {2z) v 1 



2(2^ + 3) 2z 8(2£ + 3)(2£ + 5) 4(2£ + 3)(^ + 1) 4(2^ - l)z 2 

z 6 (3£ 3 - 27^ 2 - 142^ - 136)ez 3 

_ 48(2£ + 3)(2j? + 5)(2j? + 7) ~ 48(^ + l){£ + 2){2£ + 3){2£ + 5) 

(£ 3 -m 2 + \7£- 4)e 2 £{£ - 1) {£ - 2) 2 e 3 ^ 



8(2^-l) 2 24(2£-l)z 3 
15£ 2 -(- 15£ - 11 2 



2(2£- l)(2l+ l)(2£ + 3) ' 



0(v 7 )j ; (A8) 

(A9) 
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where z — ujr and e — 2Mui, and z ~ v and e ~ v in the post New tonian expansion. The solution 4> c v can be 



obtained by changing £—»—£— 1 for £ ^ 0, 1. Incidentally, Eq. (A.8) shows the breakdown of the post-Newtonian 
expansion of i?^ -1 for ^ = and 1 at the second and third PN orders, respectively. The Wronskian of <fr" and <\>~^ v ~ 1 
becomes 

uWWtfA" l ) = — 8^WTT)WT3P e (A10) 

It should be noted that 0^ contains a fractional power of z; z v . Such a factor could be an obstacle against the power 



expansion. However, as seen from Eq. (A.4), since z v is associated with z'~ v (or z~ v with z' u ), the fractional 
power appears always in the form (z/z') v = (rjr') v (or [z! / z) u = Or' /r) v ). Therefore it may be expanded in powers 
of Or — r'), hence of r. 

To determine up to which order of (r — r') we have to expand for the 3PN calculation, we note the fact that the 
0(e 2 ) (= 0(i; 6 )) correction in v will give the factor 



Vol 2 + 151 - 11 2 
2(2e- 1)(2£ + l)(2£ + 3) 6 



in front of the powers of Or — r'). Since u> = mv t /v^ under the stationary phase approximation (which is exact for a 
circular orbit), the above term behaves as ss £ after integrating over uj and taking the m-sum. The divergence after 
summing over £ will therefore be oc (r — r')~ 2 . Thus the expansion to O (Or — r') 3 ) will be sufficient, 



15^+15^—11 
V \ ~ 2(2f-l)(2£ + l)(2£ + 3) ' 

17 



m - u 1 (^) 2 + K^) 3+o((r - /,1 »[- (Au) 



e 



2(2l-l)(2f+l)(2£ + 3) r' 2 



The coefficients of the power expansion for the radial Green function g^™^ are summarized in Appendix [TJ. We only 
note here that the leading term of the post- Newtonian expansion in gfl^ behaves as 1/r or 1/r'. 

2. Radiative part for I > 2 

To the 3PN order, the post-Newtonian orders of {j3 u , 7^} shown in Eq. ( |A6| ) imply that only the term proportional 
to 4>c( r ) ( f ) c( r ') contributes to the radiative Green function fl^Jj, 



(A12) 



The leading order term is then 



g^(ry) = » ( } 2 , + 1 + , (A13) 

which is O0v 2l+1 ) relative to the leading term in g^^ ■ Hence only the mode I — 2 is necessary to the 3PN order. 
It may be noted that under the expansion to a finite post-Newtonian order, only a finite number of £ contributes to 

(rad) 

the radiative Green function. Note also that because g\ m J is pure imaginary the radiative Green function does not 
contribute to the r-component of the reaction force in the present case of circular orbit around a black hole, whereas 
the t- and ^-components of the force are due to this part. This is in accordance with the fact that the radiative Green 
function describes the energy and angular momentum loss. 

3. I = and 1 modes 

For the modes £ = 0, 1, since the homogeneous function 4>~^ v ~ x cannot be used, we explicitly perform the post- 
Newtonian expansion of anc ^ ^up^ as done in ^,0,0. In the limit e — > 0, we have 



11 



z 2 z 4 



(0) C (0) W=ioW = l-^ + 



5040 


7x 
J ) 


z 4 z 6 
280 15120 




(1( 


z 2 i z 3 


z 4 


-y-- + 


24 ' 


!Z 3 Z 4 


i z 5 



24 120 720 



30 144 



(A14) 
(A15) 

(A16) 

(A17) 



where j n and ^ are the spherical Bessel function and the spherical Hankcl function of the 1st kind, respectively. 
Starting from the above, the expansion to the 3PN order, that is, to 0(e 2 ), is easily done to give 



A"(0)/ 



C (1) (*) = 



1 - 


z" z" 




y + 120 


l 




'l ^ 2 

^To + 


z 4 


§( 


280 


Z 


I 1 + iz- 


z 2 
2 
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5z 17z 3 
~6 + 180 



1 e 7ze 151 z 3 e 

2 I ~ "20" + 5040 



e 2 ln(z) + 0(v 7 ), 



-^ £ 2 ln(z) + 0(, 7 )) , 



iz 3 z 4 iz 5 z 6 
_ 6 _+ 24 + T20 _ 720 + 2z 



1 e 34 it 7ez 3iez 
_ _ _ _ _ + __ 



2 z 3 e 
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(A18) 
(A19) 



<t$H*) = 



—1 



e 2 

3z 2 


17 ie 2 e 3 

15z + 4z 3 


11 . 


z 2 


iz 3 z 4 


i z 5 


+ y- 




-30 + 



z" e lez 
144 + ~z 6~ 



ez + 



9e 2 



4e 3 
10 z 2 + 5I 3 



(A20) 



+ | e 2 ln(z)-^ln(z) + 0(, 7 ) 



(A21) 



The Wronskians Wi muJ — Wi muJ (<t>i^ , <jhip^ ) f° r ^ = an d 1 become 

^VK 0mw = -^-(-180 + 408«e + 935e 2 ) + O(i; 7 ), 
180 



u>W lmu = — (-225 + 311e 2 ) + 0(i; 7 ). 
225 



From the above formulas, we find the full radial Green function for £ = and 1 as 
6M 3 + 3r 3 + 3r 2 M + 4r<M 2 



9omuj{r, r') 



3r 4 



(A22) 
(A23) 



(300 M 2 r> r< + 90 r 4 + 30 r 2 r< M - 3740 M 2 r 2 + 921 r 3 M + 1320 M 2 r| ln(W>) 



180 r 3 V 

+300 Mr > r 2 +30 r 2 r 2 + 720 r 3 M ln(wr<) + 40 r 2 M 2 - 1320 M 2 ln(wr<) r 2 )w 
(r 2 +r 2 + 10 Mr < + 10 Mr>) u 3 
+ 36 Q r 2 (30 r 2 r 3 + 307 r 2 r 2 M + 15 r| + 240r 2 r 2 Mln(wr<) 



1 



+300 Mr > r 3 + 240 r 4 Mln(W<) + 3r 4 M + 68 Mr 3 r< + 212 r 4 M + 3 r 4 r<)w 4 

1 35 r 2 r 4 +7 r^ +21 r 4 r 2 



gimu>(r,r') = 



+ i (3 r< + 3 ri + 10 r 2 r 2 ) 
360 v < > > <; 5040 r< 

-10 Mr> r 2 - 18 M 2 r> r< + 18 M 3 r< - 5 r> r 3 + 5 r 3 M - 32 r> M 3 + 10 M 2 r 2 

15r| 

^tt-T (630 r 2 rM 2 - 225 r> r 4 + 1 140 r> M 2 r 2 ln(wr>) + 225 r 4 M + 900 M 2 r 3 



(A24) 
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+162 rlM 2 - 2488 M 2 r> r\ + 45 r^r 2 , - 1140 r> M 2 ln(wr<)r^ + 90 Mr^r 
+315 r 2 r 2 M - 900 Mr>r^w 2 - | (-r> r< + Mr< + Mr>)w 3 

^-3- (315 Mr 5 K - 882 Mr 2 r% - 126 r 3 r 3 + 9 r 5 >r< + 18 r5,M - 3360 r> Mln(W<) r 
+151 Mr^r< - 315 r> r^. - 504 r 3 r 2 m)u 4 - ^ir>r< (r< +r 2 ) lu 5 



7560 , < 



1 r> (105 + 189 r 2 r 4 + 27 r 4 r 2 - rg.) 



45360 r 2 



90 

w e , (A25) 



where r> = max{r, r'} and r< = min{r, r'}. 



APPENDIX B: SUMMARY OF POWER EXPANSION COEFFICIENTS 

(sym) 



In this appendix, we summarize the power expansion of the time-symmetric part of the radial Green function 
It is expanded as 



/ n-\-& ln+l \ 

^WVoHEf^'-r)^^^ , (Bl) 

where the coefficients necessary to the 3PN order calculation become 

in( _ 3) , _ (£ 3 -H£ 2 + 8£-4)£ M* 
9 tmu V,r )- 3{2 £ - l)(2£ + 1) r '3 ' ^ 

i„(- 2)(/ , = (^-l) 2 ^ ^ 2 _ (*-l) a * Mf)2 _ ^ 3 -^ 2 -^+l) M3 
yw 1 ' 0; (2£- 1)(2£+ 1) r' 2 2(2£- 1) 2 (2^+ 1) V ; (2^ - 1)(2£+ 1) r' 3 ' V ' 

in(-l), / v ^ M £ . ,., ^+1)M 2 £ . 

' (r ' ro) = 27TT- + 2(2£-l)(2^+l) ( " r )( " M) + + 8(2£-3)(2f-l)(2^+l) ( " r ) W 

, ^ 2 + 7£-4 >M)2 + ^ 3 + 5^ + 81 + 4)^ (B4) 



2(2^-l)(2£+l) v y (2£+l)(2£ + 3) r' 3 

(+0), / x 1 1 , r,2 t+l M 1 , 4 

(r ,r Q ) = -—— - -— — - — - — — — — — — — — - — (ur ) 



2£+l 2(2£ - 1)(2£ + 1) v ' 2£ + 1 r> 8(2£ - 3){2£ - l)(2£ + 1) 



21(21- l)(2£ + l) v M ; (2£+ l)(2£ + 3) r' 2 48(2£ - 5)(2£ - 2>){2l - l)(2l + 1) 

■K) 3 (wM) 



3£ 3 • 36/ 2 7!)' • 2 1 



24(€ - l)£{2£ - 3)(2£ - \){2l + 1) 

■(wJlf) 



1488^ 6 + 7104£ 5 + 10608^ 4 + 1940£ 3 - 6825^ 2 - 2321£ + 1740 , _ n2 



6(2*- 1) 2 (2£+ 1) 3 (2^ + 3) 2 

2(15£ 2 + 15£-11) . -„ 2 , . , ^ 4 + 9^ 3 + 29^ 2 + 39<? + 18 M 3 

+ (2^-l)(2£+l)2(2^ + 3) ( " M) ln(r ° /r } 3(21 +l)(2l + 3) 7^ ' (B5) 

inf+i)/ / x ^ 2 -5£-10 , , w w , £ 2 -5^-10 , ,, 3/ 

^ (r ' ro) = ' 2(^+l)(2£+l)(2£ + 3) ( " r )(WM) " 4(£+l)(2£-l)(2^+l)(2^ + 3) ( " r } ( " M) 

6(15£ 2 + 15£ - 11) r' , 2 ^ 2 -5^-10 . , r , 2 

■(wM) 2 - — — — — — — — (u>M) 2 , (B6) 



(2^-l)(2£+l) 2 (2^ + 3)r v ' 2(2£+ l)(2£ + 3) 
(+2) (r',r ) = — l ———(ur'f + — .. / — - - (o>/) 4 + — — — - (wr 7 ) (uM) 



2(2£+l)(2^ + 3) v ' A{2£- l){2£ + l)(2£ + 3) v y 2(2^ + l){2£ + 3) 



16(2^ -3)(2^- l){2£+ l)(2£ + 3) v ' U{2£ - l)(2£ + 1)(2£ + 3) 
3(15l 2 + 15l-ll) r 72 (f+l)(£ + 2) 2 2 

+ (2^-l)(2£+ 1)2(2^ + 3) rf M/) " 2(2£ + 1)(2^ + 3) 2 ( " M) ' (B?) 
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+ 8(2l-3)(2f-l)(2^ + l) ^^' ( B19 ) 
7^~V,n>) = 0, (B20) 
iZL V^o) = ~ 48(2 ^ _ - 1)(2 ^ _ 1)(2£ _ 3) ( w/ ) 6 ■ (B21) 



APPENDIX C: GENERATING FUNCTION OF M SUM 

In thi s app endix, we describe a method to sum over the m-modes of spherical harmonics for arbitrary i that appears 
in Eq. ( |3.1S| ). Specifically, the m-sum we need to evaluate takes the form, 

i 

£ m w |y £m (7r/2,0)| 2 , (CI) 

m=-l 

where N is a non-negative integer. To perform the above summation, we introduce the generating function, 

i 

T i {z)= J2 e mz \Y em (n/2,0)\ 2 . (C2) 



Then the sum (CI) may be evaluated as lim d^Ti(z). 

z— >0 



First note that Y^ m (7r/2, 0) is non-zero only when I — m is an even integer. So putting 2n = I — to, we have 



Tf(z) _y V i-2n), 2'+l (2rQ! / (2l-2n-l)!! \ 
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E 



4tt 4 £ ^ r(n + l)r(^-n+l) 2 n\ 

2l+ir(l + l/2)r(-l+l/2) tz T(n + l/2)r(n - I) e~ 2 " 2 
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+-;e-^ , (C3) 



4^ 2 T(£+l) z 1 V2' ' 2 

where 2-F1 is the hypergeometric function, and we have temporarily performed the analytic continuation of I to a 
complex number from the second line to the third line. 

Using the relation between the hypergeometric functions as 

(^4*)='--> , r(<+ ^r^ -^H' <«> 



2 Fi (j, -4 2 ; x 

for a non-negative integer I, we have 



This can be easily expanded to arbitrary order of z. For example, to 0(z 6 ) we have 

2£+ 1 / / £(g+m 1 a ^(£+l)(3l 2 + 3l-2) ^ 1 4 



^+l)(5J? 4 + 10f 3 -5£ 2 -l(W + 8)\ 1 (1 



16 / 6! 



-z b + 0{z») . (C6) 
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